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Using  the  complex  Toda  chain  (CTC)  as  a  model  for  the  propagation  of  the  TV'-soliton  pulse  trains 
of  the  nonlinear  Schrodinger  (NLS)  equation,  we  can  predict  the  stability  and  the  asymptotic  behav¬ 
ior  of  these  trains.  We  show  that  the  following  asymptotic  regimes  are  stable:  (i)  asymptotically  free 
propagation  of  all  N  solitons;  (ii)  bound  state  regime  where  the  N  solitons  move  quasi-equidistantly; 
and  (iii)  varioxis  different  combinations  of  (i)  and  (ii).  On  the  example  of  TV  =  3  we  show  how  the 
CTC  model  can  be  used  to  determine  analytically  the  set  of  initial  soliton  parameters  corresponding 
to  regime  (ii).  We  compare  these  analytical  results  against  the  corresponding  numerical  solutions 
of  the  NLS  and  find  excellent  agreement  in  most  cases.  We  concentrate  on  the  quasi-equidistant 
propagation  of  all  N  solitons  because  it  is  of  importance  for  optical  fiber  soliton  communication. 
We  check  numerically  that  such  propagation  takes  place  for  TV  =  2  to  8.  Finally  we  propose  realistic 
configurations  for  the  sets  of  the  amplitudes,  for  which  the  trains  show  quasi-equidistant  behavior 
to  very  large  run  lengths. 


Submitted  to  Phys.  Rev.  Lett. 


I.  INTRODUCTION 

One  of  the  important  problems  in  optical  fiber  soliton 
communication  is  to  achieve  as  high  of  a  bit  rate  as  pos^ 
sible.  In  order  to  do  this,  one  needs  to  be  able  to  pack 
the  solitons  into  as  short  of  a  space  as  possible.  However, 
if  the  solitons  are  too  close  together,  then  their  mutual 
(linear  and  nonlinear)  interactions  can  cause  them  to  col¬ 
lide  and/or  separate,  thereby  corrupting  the  signal.  The 
current  solution  of  this  problem  is  simply  to  require  each 
soliton  to  be  sufficiently  far  apart  from  all  others  (usu¬ 
ally  6  or  so  soliton  widths)  so  that  such  interactions  can 
be  totally  neglected.  However,  at  the  same  time,  it  was 
predicted  [1,2]  and  experimentally  confirmed  [3]  that  for 
certain  values  of  relative  soliton  parameters,  this  separa¬ 
tion  can  be  reduced,  and  at  the  same  time,  still  maintain 
signal  integrity.  Our  main  purpose  here  is  to  analytically 
and  numerically  detail  the  soliton  parameter  regime,  in¬ 
side  of  which,  signal  integrity  can  be  maintained.  In  par¬ 
ticular,  we  are  interested  to  determine  how  one  may  use 
this  inter-soliton  interaction  for  stabilizing  a  soliton  train. 

Any  communication  signal  will  be  composed  of  ’’ran¬ 
dom”  combinations  of  O’s  and  I’s.  It  cein  also  be  viewed 
as  being  composed  of  a  random  collection  of  TV-soliton 
trains,  with  varying  widths  of  O’s  between  them.  Thus  it 
is  then  adequate  for  us  to  simply  analyze  the  stability  of 
individual  A-soliton  trains,  for  finite  N.  This  we  will  do 
and  will  study  a  nonperiodic  and  finite  train  (chain)  of 
soliton  pulses  by  both  analytical  and  numerical  methods 
using  and  developing  the  ideas  in  [4,5]. 

The  basic  model  and  description  of  AT-soliton  trains  in 
optical  fibers  is  provided  by  the  nonlinear  Schrodinger 
(NLS)  equation  and  its  perturbed  version: 

iut  +  +  \u\^u{x,  t)  =  iR[u],  (1) 


where  0k,o{^)  ~  ~  ^k,o)  ^^id  6^,0^  sfc.o- 

and  2/ifc,o  are  the  initial  phase,  position,  amplitude  and 
velocity  respectively  of  the  /c-th  soliton;  by  2vo  and  2fio 
we  mean  the  averages  of  the  initial  amplitudes  and  ve¬ 
locities. 

An  effective  method  for  studying  the  interaction  of 
such  trains  of  soliton  pulses  was  first  proposed  by  Karp- 
man  and  Solov’ev  (KS),  for  the  simplest  nontrivial  case 
of  a  2-soliton  interaction  [12].  For  further  developments 
and  analysis  for  different  physically  important  perturba¬ 
tions  can  be  found  in  [13,14]  and  the  references  therein; 
for  a  review  see  Ref.  [7].  The  KS  method  is  based  on 
the  adiabatic  approximation.  It  is  valid  for  any  collec¬ 
tion  of  well  separated  solitons,  such  that  their  mutual 
interactions  will  lead  to  a  slow  deformation  in  the  soliton 
parameters.  As  a  result  the  soliton  interactions  should 
be  described  by  a  dynamical  system  for  Mt  ^tnd  6(. 

Gorshkov  [15]  and  Arnold  [16],  have  conjectured  that 
an  infinite  train  of  out-of-phase  soliton  pulses,  with  equal 
amplitudes  and  velocities  could  be  described  by  the  real 
Toda  chain: 


^  f,qk^\-qk  _  pqk-qk~\ 

dT^ 


(3) 


where  r  =  4uot,  A:  =  0,  ±1,  ±2, . . .  and  qk  are  real  func¬ 
tions  related  to  the  soliton  positions.  This  system  will 
be  referred  to  as  the  real  Toda  chain  (F{TC).  The  valid¬ 
ity  of  this  description  for  finite  number  of  solitons  (see 
also  below)  has  been  verified  also  numerically  [17]. 

Recently  in  Refs.  [4,5,17],  the  Karpman-Solov’ev 
method  was  extended  to  iV-soliton  pulses,  and  then  with 
additional  approximations,  was  reduced  to  the  complex 
Toda  chain  equations  (CTC)  (3)  with  N  sites.  The  cor¬ 
responding  system  of  equations  is  (3),  but  with  k  = 
1,...,^  and  =  o.  The  complex  valued 

functions  qk{t)  are  related  to  the  soliton  parameters  by: 


This  equation  describes  a  variety  of  wave  interactions, 
including  solitons  in  nonlinear  fiber  optics  [1,6-11]  and 
spatial  solitons  in  nonlinear  refractive  media  [9]. 

The  inverse  scattering  method  [6]  allows  one  to  solve 
exactly  Eq.  (1)  when  R[u]  =  0  and  to  calculate  explicitly 
its  iV-soliton  solutions.  However  for  our  purposes,  this 
method  is  impractical  for  two  reasons.  First,  there  are 
important  perturbations  of  this  system  which  have  no 
explicit  solutions.  Second,  an  approximate  method  can 
serve  much  better  than  an  exact  approach  since  the  AT- 
soliton  trains  here  are  rather  special  and  can  only  be 
approximated  by  iV-soliton  solutions.  Such  trains  are 
actually  sums  of  1-soliton  pulses,  which  are  spaced  almost 
equally  apart,  have  almost  equal  amplitudes,  and  move 
with  essentially  the  same  velocity.  More  specifically,  they 
are  the  solutions  to  Eq.  (1)  with  R[u]  =  0  satisfying  the 
following  initial  conditions: 

~^cosh(2i/fc(a:-eM)’ 


Qk+I  -Qk  =  -2i'o(Cfc+i  -  ^fc)  +  ln4i/o 

-h  i  (tt  -f  2^(^fc+i  -  Cfc)  “  (4+1  “  4)) .  (4) 

Thus  the  problem  of  determining  the  evolution  of  an  NLS 
AT-soliton  train  has  been  reduced  to  the  problem  of  solv¬ 
ing  the  CTC  for  N  sites.  CTC,  like  RTC  is  integrable 
and  we  may  use  the  special  techniques  valid  for  integrable 
lattices,  see  [18-20]. 

Our  main  results  are  the  following: 

(i)  we  show  by  analyzing  the  exact  analytic  solutions  of 
CTC,  that  it  has  several  qualitatively  different  classes  of 
asymptotic  regimes.  Besides  the  asymptotically  free  mo¬ 
tion  (which  is  the  only  one  possible  for  the  RTC  [19,20]), 
CTC  allows  also  for:  (a)  bound  state  regime  when  all 
the  N  particles  move  quasi-equidistantly;  (b)  all  possi¬ 
ble  intermediate  regimes  when  one  (or  several)  group(s) 
of  particles  form  bound  state(s)  and  the  rest  of  them 
go  into  free  motion  asymptotics.  In  addition  to  these 
relatively  stable  regimes  of  motion  there  are  also  less  sta¬ 
ble  regions  in  the  space  of  soliton  parameters,  where  one 
regime  switches  into  another  one.  There  one  can  find 
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(c)  singular  solutions,  and  (d)  various  types  of  degener¬ 
ate  solutions. 

(ii)  we  compare  the  predictions  from  the  CTC  model 
with  the  numerical  solutions  of  the  NLS,  that  bound 
states  regimes  indeed  take  place  in  the  soliton  interac¬ 
tions  and  are  very  well  described  by  the  CTC  model. 
Our  analytic  approach  allows  us  to  predict  the  set  of 
initial  parameters,  for  which  these  asymptotic  regimes 
takes  place.  We  put  special  stress  on  the  bound  state 
and  quasi-equidistant  regimes  since  such  solutions  would 
be  the  desirable  in  long  range  fiber  optics  communica¬ 
tions. 


n.  ASYMPTOTIC  REGIMES  OF  THE  CTC 

As  in  [4,5],  one  can  generalize  the  RTC  [18-20]  to  the 
complex  CTC  case.  We  list  the  four  most  important 
points  concerning  this  below: 

51)  The  CTC  Lax  representation  is  the  same  as  for  the 
RTC:  L  =  [B,  Lj,  where 

N 

L  =  ^  {HEkk  +  <^kiEk,k^i  H-  Ek^i^k)) »  (5) 

fc=i 

Here  ak  =  ie(9fc+i-9k)/2  ^nd  bk  =  The 

matrices  {Ekn)pq  =  ^kp^nqi  and  {Ekn)pq  =  0  whenever  p 
or  q  becomes  0  or  AT  +  1. 

52)  The  integrals  of  motion  in  involution  are  provided 
by  the  eigenvalues,  of  Lo  =  i>(T  =  0). 

53)  The  solutions  of  both  the  CTC  and  the  RTC  are 

determined  by  the  scattering  data  5lo  of  Lq-  If  Cfc  #  Cj 
for  k  ^  j,  then  Here  are  the  first 

components  of  the  corresponding  eigenvectors  of  Lq? 
which  are  fixed  (up  to  an  overall  sign)  by  the  condition 

=  I?  [see  [19,20]].  S'x.Q  uniquely  determines 

both  Lq  and  the  solution  of  the  CTC. 

54)  Lastly,  Ck  uniquely  determine  the  asymptotic  be¬ 
havior  of  the  solutions  of  the  CTC  and  can  be  calculated 
from  the  initial  conditions.  We  will  use  this  fact  to  de¬ 
scribe  the  different  classes  of  asymptotic  behavior. 

In  addition  to  the  dynamical  variables  becoming  com¬ 
plex  valued,  there  are  other,  important  differences  be¬ 
tween  the  RTC  and  CTC,  and  their  properties.  For  the 
RTC,  one  has  that  [19,20],  both  the  eigenvalues,  Ckj  and 
the  coefficients,  r*.,  are  always  real-valued.  Moreover, 
one  can  prove  that  Ck  ^  Cj  for  j-  As  a  consequence 
the  only  possible  asymptotic  behavior  in  the  RTC  is  an 
asymptotically  separating,  free  motion  of  the  solitons. 

For  CTC  not  only  the  dynamical  variables  qk ,  but  also 
Cfc  =  /Cfc  4-  irjk  and  rk  take  complex  values.  The  collec¬ 
tion  of  eigenvalues,  still  determines  the  asymptotic 
behavior  of  the  solitons.  In  particular,  it  is  Kk  that  de¬ 
termines  the  asymptotic  velocity  of  the  A:-th  soliton.  For 
simplicity  and  without  loss  of  generality  we  assume  that: 
trLo  =  0;  Cfc  ^  Cj  for  k  ^  j  (this  does  not  necessarily 


mean  that  Kk  ^  Kj)\  and  that  the  Kk's  are  ordered  as: 

^  ^2  £  ■  •  •  £  Then  in  any  train  of  solitons,  there 
are  three  possible  general  configurations: 

Dl)  Kk  ^  Kj  for  k  ^  i.e.  the  asymptotic  velocities 
are  all  different.  Then  we  have  asymptotically  separating, 
free  solitons,  see  [4,5] 

D2)  Ki  —  K2  =  . .  •  =  i-e.  all  N  solitons  wiW  move 
with  the  same  mean  asymptotic  velocity,  and  therefore 
will  form  a  “bound  state” .  The  key  question  now  will  be 
the  nature  of  the  internal  motions  in  such  a  bound  state. 
In  particular,  one  would  want  any  two  adjacent  solitons 
to  move  quasi-equidistantly. 

D3)  One  may  have  also  a  variety  of  intermediate  situa¬ 
tions  when  only  one  group  (or  several  groups)  of  particles 
move  with  the  same  mean  asymptotic  velocity;  then  they 
would  form  one  (or  several)  bound  state(s)  and  the  rest 
of  the  particles  will  have  free  asymptotic  motion. 

Obviously  the  cases  D2)  and  D3)  have  no  analogies  in 
the  RTC  and  physically  are  qualitatively  different  from 
Dl).  The  same  is  also  true  for  the  special  degenerate 
cases,  where  two  or  more  of  the  life’s  may  become  equal. 
These  cases  will  be  considered  elsewhere.  Another  type 
of  solutions  of  the  CTC  which  should  be  dealt  with  sep¬ 
arately  are  the  singular  solutions,  see  e.g.  [22]. 

We  will  say  that  the  motion  is  “quasi-equidistant”  if 
Afc  1,  where 

Ak  =  (max(^fc+i  -  ffc)  -  -  ik))  Ao,  (6) 

then  the  solitons  will  not  asymptotically  separate,  but 
instead  will  slowly  oscillate  with  some  small  amplitude, 

iFrom  here  on,  we  will  denote  the  various  sets  of  ini¬ 
tial  conditions  for  the  AT-soliton  trains  by  quadruplets 
Ar|ro|A|Ph,  where  N  is  the  number  of  pulses,  ro  is  the 
distance  between  neighboring  pulses  and  A  and  Ph  stand 
for  the  configuration  of  initial  amplitudes  and  phases  re¬ 
spectively.  In  our  runs,  listed  in  the  tables,  the  pulses  are 
initially  equidistant  (i.e.  ""  CkyO  =  ro  =  8,  which 

corresponds  to  4  pulsewidths),  and  the  initial  velocities, 
2/Xfco  are  all  vanishing.  We  use  here  several  configurations 
of  amplitudes: 

Ai:  the  amplitudes  are  monotonic,  i.e.  Az/q  = 
z/fc,o  =  2di  and  such  that  the  average  amplitude  equals 
2i/o;  this  we  denote  by  the  number  2Az/o  lO^; 

A2:  configurations  with  two  different  alternating  am¬ 
plitudes  z/1,0,  1^,0,  1^1,0,  •  •  • ; 

A3(4):  configurations  with  three  (or  four)  different  al¬ 
ternating  amplitudes  z/i,o>  •  •  • ,  *^3(4);0)  *^1,0  •  •  • ;  these  we 
denote  by  “3sw”  (or  “4sw”)  assuming  that  is  the 
smallest  one. 


m.  THE  SOLITON  BOUND  STATES 
A.  Analytical  results 

For  the  sake  of  brevity  we  present  here  some  of  the 
analytical  results  only  for  iV  =  3;  it  is  possible  to  extend 
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these  results  to  any  A^ 

We  assume  that  ~  center  of 

mass  is  fixed  at  the  origin.  This  is  compatible  with 
tr  Lo  =  0.  The  initial  soliton  parameters  determine 
Lo  through  (4)  and  (5)  and  the  regimes  of  propagation 
are  determined  by  the  eigenvalues  of  Lo-  In  particular, 
the  three  solitons  form  a  bound  state  if  the  eigenvalues 

=  irjf^  are  purely  imaginary. 

The  explicit  analytic  solution  to  the  N  =  S  CTC  is 
well  known  [19,20].  Choosing  the  ’’symmetric”  one  for 
which  <72  =  0,  =  -<73  and  =  (2  -  k)iT]i  we  find  that 

it  is  periodic  with  period  —  7r/{2i/o7?f )  where 

2rit  =  Y/(Ai/o)2±(Ai/cr,3)''  (7) 

Az/cr,3  =  2y/2eo  and  eo  =  The  choices  for  the 

soliton  parameters  here  are  Ai  and  Ph  =  {0, 62,07  0}» 
and  the  signs  plus  and  minus  correspond  to  ^2,0  =  0  and 
^2  0  =  TT  respectively.  For  this  solution  we  find 

where  zq  =  |Az/o|/Ai/cr,3-  Note  that  similar  results  for 
Af  hold  also  for  the  N  —  2  case  with  the  only  difference 
being  that  Az/cr,3  should  be  replaced  by  4eo.  The  motion 
will  be  quasi-equidistant  if  Af  <  1.  The  formulas  (8) 
show  that  the  increase  of  |Ai^o|  diminishes  Af .  Another 
way  to  diminish  A\  for  fixed  Ai/q  and  ro,  is  to  increase 
the  average  amplitude  i/q* 

The  singular  behavior  of  Af  for  Az/q  0  corresponds 
to  a  singular  solution  of  the  CTC;  the  numeric  solution 
of  the  NLS  shows  that  the  solitons  do  not  collide,  but  do 
come  rather  close  to  each  other  at  the  values  of  t  where 
the  CTC  solution  develop  singularities.  We  will  come 
back  to  this  question  in  Subsection  IILB  below. 

On  the  other  hand  the  singularity  of  A^  for  Az/q  = 
Ai/cr,3  corresponds  to  the  fact  that  at  this  critical  point 
the  quasi-equidistant  regime  switches  over  into  the  free 
motion  regime.  Indeed,  as  we  will  see  below,  for  Az/q  < 
Az/cr  3  we  have  regime  (i)  and  the  distance  between  the 
solitons  grows  infhiitely,  while  for  Az/q  >  Ai/cr,3  we  get 
regime  (ii)  and  a  possible  quasi-equidistant  behavior. 

Let  us  now  describe  the  particular  choices  of  the  soliton 
parameters,  which  lead  to  the  regimes  described  above. 
Such  analysis  must  be  based  on  the  solution  of  the  char¬ 
acteristic  equation  for  Lq,  which  for  Af  =  3  with  tr  Lo  =  0 
and  generic  choice  for  A  and  Ph  is: 


C“  +  Cp  +  9  =  o, 

(9a) 

—  (dj  +  dj  +  <^)  —  a?  —  02* 

(9b) 

z  z 

-  [a^ds  +  a^di)  +  —^1^2^31 

(9c) 

where  dk  =  ^  and 

Now  we  can  use  Cardano  formulas  to  evaluate  the  roots 
ijfc  which  determine  the  asymptotic  regime  predicted  by 
CTC.  For  brevity  and  simplicity  we  limit  ourselves  to 
two  special  choices  of  the  initial  amplitudes.  Skipping 
the  details,  we  find  that  the  bound  state  regime  occurs 
for  each  of  the  following  choices  of  the  soliton  parameters. 

Case  Ai. 

ii.a)  Ph  =  {0,0,0}  and  |Az/o|  >  0: 

ii.b)  Ph  =  {0,  drTT, 0}  and  |Az/o|  >  2y/2eo  =  Ai/cr.3; 

ii.c)  Ph  =  {0,<5o,<5o}  and  Ph  =  {0,0,  <5o}  where  <5o  = 
±7r  and  |Az/o|  >  2  • 

Case  A2. 

ii.d)  Ph  =  {0,0,0}  and  |Az/o|  >  0; 

ii.e)  Ph  =  {0,  d::7r,0}  and  [Az^ol  >  4\/2eo; 

ii.f)  Ph  =  {0,52,0^0}  and  \Ai/o\  >  0  if  00362,0  >  0; 
if  cos  62,0  <  0  then  the  bound  state  regime  holds  for 
|Ai/o|  >  4yf— 2cos  62,0^0; 

As  one  may  expect  the  asymptotic  behavior  depends 
very  much  on  the  initial  choice  of  phases.  Thus  in  the 
cases  ii.a)  and  ii.d)  the  bound  state  regime  is  obtained 
for  any  value  of  Az/q,  while  in  all  other  cases  this  regime 
is  entered  only  if  |Az/o|  is  larger  than  some  critical  value, 
which  of  course,  depends  on  the  initial  parameters,  com¬ 
pare  e.g.  ii.b),  i.c),  ii.e)  and  ii.f). 

B.  Comparison  Between  CTC  and  NLS 

In  this  subsection  we  compare  the  analytical  results 
obtained  above  with  the  numerical  solutions  of  the  un¬ 
perturbed  NLS  (ii[u]  =  0)  (1)  with  initial  conditions  (2). 

In  tables  I,  II  and  III  we  list  the  values  of  maxAfc 
obtained  from  the  niimerical  solution  of  the  NLS  equation 
with  the  corresponding  CTC  solution.  It  is  seen  that  the 
agreement  is  very  good  and  improves  with  the  increasing 
of  Az/q  and  Tq,  In  addition  we  show  that  maxAfc  (6) 
indeed  has  very  low  values,  which  is  a  reason  to  call  such 
regime  quasi-equidistant. 

We  find  an  excellent  match  between  CTC  and  NLS  for 
the  description  of  all  soliton  parameters  in  most  cases. 
The  exceptions  in  this  respect  happen  for  those  sets  of 
soliton  parameters  at  which:  a)  CTC  has  singular  or  de¬ 
generate  solutions  and  b)  transition  from  one  regime  to 
another  takes  place.  The  choice  of  Ai  with  Az/q  =  0  and 
Sfc  Q  =  0  leads  to  singular  solution  of  CTC  for  any  AT.  On 
the  other  hand  the  NLS  solutions  are  always  analytic  and 
never  have  singularities.  Our  numerical  checks  show  that 
such  sets  of  initial  conditions  correspond  either  to  solitons 
collisions  or  soliton  coalescence  depending  on  the  value 
of  ro.  In  the  regions  where  this  happens,  the  adiabatic 
approximation  is  no  longer  valid  and  the  CTC  does  not 
match  with  the  numeric  solutions.  For  AT  =  3  and  ro  =  8 
the  first  coalescence  takes  place  at  ti  (T^/2)|ai/o=o 
and  the  next  ones  tend  to  repeat  periodically  with  pe¬ 
riod  very  close  to  r3'lAi/o=o  (7).  The  quasi-equidistant 
propagation  of  the  solitons  is  maintained  for  t  <  0.9ti . 
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For  2Ai/o  >  0.1  we  see  that  in  the  quasi-equidistant 
regime  the  soliton  positions  display  periodic  behavior 
with  small  periods  (much  less  than  ti)  and  with  small 
amplitudes  of  oscillations  given  by  In  several  occa> 
sions  we  extended  the  runs  to  lengths  of  t  =  600  and  996 
getting  the  same  type  of  behavior  and  the  same  results 
for  Ak-  Thus  we  conclude  that  these  sets  of  initial  pa- 
rameters  lead  to  a  very  stable  behavior.  In  particular, 
for  Ai  configurations  with  2Ai^o  =0.2  and  tq  =  8  we 
have  an  equidistant  behavior  with  an  error  ^  10%. 

The  bound  state  regime  is  characteristic  for  Ai  trains 
with  any  N.  However  the  Ai  configurations  become  im¬ 
practical  if  one  needs  a  very  long  sequences  of  bytes  all 
equal  to  one.  so  we  should  explore  other  possibilities. 

One  of  them  is  to  use  a  soliton  train  of  in-phase 
solitons,  with  A2  amplitude  configurations  [2],  where 
2i/i  =  2i/3  =  1.0,  2i/2  =  2z/4  =  1.25  and  the  rtm  length 
is  about  140.  This  idea  was  verified  experimentally  by 
[3J  where  20  Gbit/s  single  channel  soliton  transmission 
over  1 1  500  km  using  alternating  amplitude  solitons  was 
reported.  To  check  and  develop  this  idea  we  did  a  se¬ 
ries  of  runs  of  N  in-phase  soliton  trains  with  A2,  A3 
and  A4  amplitude  configurations  with  A  =  3  to  8,  see 
Tables  II  and  III.  For  the  A2  and  A3  configurations  we 
find  that  the  quasi-equidistant  propagation  takes  place 
for  t  <  0.9Tqed  with  ~  250.  For  larger  values  of 

t  some  of  the  solitons  come  rather  close  to  each  other. 
CTC  also  predicts  the  same  type  of  behavior  but  with 
larger  value  for  ~  460.  This  is  the  reason  why  in 

Table  II  the  numerical  data  from  NLS  shows  non  equidis¬ 
tant  propagation  when  considered  for  run  length  300.  We 
extended  some  of  these  runs  to  lengths  of  600  and  996. 
The  results  for  iV  =  3  show  a  structure  close  to  a  peri¬ 
odic  one  with  a  period  of  about  2T^^.  So  if  the  goal 
is  to  achieve  a  quasi-equidistant  propagation  to  lengths 
t  <  then  such  configurations  can  be  used. 

For  larger  run  lengths  other  configurations  must  be 
used,  see  also  [1].  The  next  possibility  which  we  explored 
is  the  A4  configurations  of  soliton  amplitudes.  The  re¬ 
sults  are  collected  in  Table  III  for  A  =  4  to  8.  We  see  a 
substantially  different  picture.  Both  the  NLS  and  CTC 
show  a  quasi-equidistant  behavior  with  the  same  small 
values  of  Ak  like  in  the  last  column  of  Table  III  to  run 
lengths  of  1 200;  for  the  CTC  this  behavior  persists  even 
to  lengths  of  12000. 

IV.  CONCLUSIONS 

A  method  for  the  description  of  the  asymptotic  behav¬ 
ior  of  the  AT-soliton  pulse  trains  of  the  NLS  equation  is 
proposed,  based  on  the  CTC  model  for  the  soliton  inter¬ 
action.  It  describes  correctly  several  qualitatively  differ¬ 
ent  classes  of  asymptotic  regimes.  Several  sets  of  soliton 
parameters  have  been  described  for  which  the  propaga¬ 
tion  is  quasi-equidistant.  Such  behavior  with  a  conve¬ 
niently  low  value  for  Ak  can  be  achieved  by:  a)  taking 


soliton  trains  with  Auq  large  enough;  b)  increasing  the 
value  21^0  of  the  average  amplitude;  and  c)  increasing  the 
distance  ro  between  the  neighboring  solitons. 

The  critical  values  of  the  soliton  parameters,  for  which 
one  regime  switches  over  to  another  one  have  been  eval¬ 
uated.  We  find  that  near  these  critical  values,  the  match 
between  CTC  and  NLS  becomes  worse,  as  expected. 

The  CTC-model  provides  one  with  a  tool  for  construct¬ 
ing  sets  of  initial  data  for  the  A'-soliton  trains  which 
will  possess  a  given  asymptotic  trait,  determined  by  the 
eigenvalues  Qk  of  Lq. 

The  monotonically  increasing  amplitudes  cannot  be 
used  in  case  one  needs  trains  with  larger  number  of  soli¬ 
tons.  That  is  why  we  investigated  A2,  A3  and  A4  ampli¬ 
tude  configurations,  for  which  the  propagation  is  quasi- 
equidistant.  We  find  that  the  A2  and  A3  amplitudes  con¬ 
figurations  may  provide  quasi-equidistant  propagation  to 
run  lengths  of  the  order  of  The  A4  configurations 

with  four  different  amplitudes  show  quasi-equidistant  be¬ 
havior  to  much  larger  run  lengths. 
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TABLE  II.  The  value  of  max  At  from  NLS  and  CTC  for 
ro  —  8  and  A2  configuration;  by  “si”  we  have  denoted  trains 
with  2i/i  =  1.0  and  21^2  =  1.25  and  “Is”  means  2i/\  =  1.25  and 
2i/2  =  1.00.  Run  lengths  equal  300. _ _ 


NLS 

CTC 

NLS 

CTC 

NLS 

CTC 

3|8|sl|0 

318|ls|0 

4|8|sl|0 

0.62 

0.08 

0.012 

0.026 

0.42 

0.078 

4|8|ls|0 

5|8|sl|0 

5|8|lsl0 

0.42 

0.078 

0.23 

0.058 

0.59 

0.031 

6|8|sl|0 

6|8|ls|0 

7|8|sl|0 

0.27 

0.039 

0.28 

0.039 

0.26 

0.073 

7|8|ls|0 

8|8|sl|0 

8|8|ls|0 

0.27 

0.035 

0.25 

0.04 

0.31 

0.04 

TABLE 

value  of  maxAjt 

III.  The 

fi'om  NLS  and  CTC  for  ro  =  8  and  A3, 

and  A4  configurations:  3sw  — »  {0.85, 1.00, 1.15,0.85, . . .}  and 
4sw  — ►  {0.85,1.00,1.15,1.30,0.85,...}.  Run  lengths  equal 

1200. 

NLS  CTC 

NLS  CTC 

NLS  CTC 

4|8|3sw|0 

5|813sw|0 

5|8|4sw|0 

0.19  0.07 

0.47  0.07 

0.06  0.04 

6|813sw|0 

6(8|4sw|0 

7|8|3sw|0 

0.46  0.20 

0.10  0.04 

0.79  0.44 

7|8|4sw|0 

8|8|3sw|0 

8|8|4sw|0 

0.09  0.04 

0.84  0.33 

0.09  0.03 

« 


TABLE  I.  The  value  of  max  Ak  from  NLS  and  CTC  for 
ro  =  8  and  Ai  configuration.  Run  lengths  equal  300. _ 


NLS 

CTC 

NLS 

CTC 

NLS 

CTC 

2|8|10|0 

21811510 

21812010 

0.055 

0.054 

0.029 

0.026 

0.019 

0.016 

31811010 

31811510 

31812010 

0.099 

0.03 

0.053 

0.014 

0.031 

0.008 

41811010 

41811510 

4|8|20|0 

0.14 

0.13 

0.085 

0.051 

0.064 

0.026 

51811010 

518|1510 

58 

12010 

0.20 

0.09 

0.13 

0.026 

0.11 

0.019 
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